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Th® p reseat work involves the seedy of thermal 
stresses in thin shallow shells* Chapter t gives the 
introduction about the t hemal stresses and msst sad 
type of work done in thin shallow shells. Chapter 11 deals 
With the variation of total potential energy in deriving 
the basic enactions for general thin shells end possible 
boundary conditions In thin shells. Chapter III deals 
with the asking of nonhomogeaoous boundary conditions to 
homogeneous conditions fey Hindi in end Goodman technique. 
This chapter gives the solutions for two parallel edges 
simply supported end two other edges supported in my 
manner. Solutions have heea given for (a) all four edges 
simply sapper tod end (b) two parallel edges simply suppo- 
rted and remaining two edges clusped, where constants 
have been obtained with the help of Hoang's Tbchaique. 
m example Illustrates the method. 

A flew chart end a computer programme for (a) ail 
edges simply supported and tb) two parallel edges simply 
supported and Sustaining two edges clasped, are gives in 
Appendix (B) and <C) respectively. Tables sad graphs are 
presented for deflections end moments in tom of dimen- 
sionless parameters for e/b ratio from 1.0 to S,0 at 
interval of O.S and for « 2 /Wi retie from £.0 to SO.O 
at interval of 6.0 for elliptic paraboloids. 
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Hendlnn end twisting moments per unit dint cm oe 

in shell. 

a£ / f1 /* fT.s* da 
-h/£ 

Oik; / h/S 2. da 

-h/2 

components of intensity of surface loading? 
'Temperature rise. 

Normal end shearing forces per unit distance 
In middle surface of shell* 

Components of displacements et distance u from 
middle plane* 

Coefficient of linear expeneion 
Curvilinear Coordinates 
formal components of stress parallel to 
@l« pg end 0-axie 

Shearing stress components in curvilinear 
coordinates. . 
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In ©ado rat design, there is a let of atiantj on 
towards the themal stresses chiefly because of many 


are also Introduced in oases where no temperature gradient 
exist# if one pert consists of several neterials having 
different coefficient of expansions, or If free contraction 
or expansion is prevented by external constraints. At an 


eeiying with depth only will bead into a served plate, Theiw 
ml stresses trill he let rodeoed in a free pitta made of too 
astsrislii having different coefficient of expansions, with 
uniform tojaperatnre distribution, sm& the real stresses are 










s 

also set in constrained edges plate having uni fern 
temperature distribution* 

i ** 

The basic so orgy e^ueti on derived frets *Le*r 

of conservation of energy* tender the action of bent is 
given as 

K r #u *Pc £ T -i <3X*£^) ^-T 0 £ kk CD 

This is a coupled eqttotiea, which royui rod simultaneous 
determination of temperature distribution sod deformations* 

This boundary value problem has considerable mathe- 
matical difficulty, since it combi oes the theories of heat 
conduction and elasticity under transient conditions. bat 
fortunately in most of engineering problems, it is possible 
to introduce certain simplifying assumptions without signi- 
ficant error* The principal simplifying assumption is omi- 
ssion of mechanical coupling term in energy equation and 
giving rise some fora of Fourier heat conduction equation, 
which is independent of stresses nod deformations. It vari- 
ations ef strain In a body are produced by some external 
agency ■ these variations of strain ore accompanied hp vari- 
ations in temperature and consequently flew ef heat. Thin 
gives rise to as increase of entropy end corresponding in- 
crease of energy stored in s mechanically irrecoverable 
manner and this is theme# las tic dissipation, which req- 
uires far its study the use ef coupled heat equation. The 
deformations dee to external lead# are accompanied by only 
small changes in temperature* It is then reasonable to 
calculate these deformations without t eking account of 
* lumbers refer to entries In the list of references. 
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thermal expansion* Similarly, if strain* era produced 
la a body fay ooe-uelfom tamps raters distri ballon, it 
Is clear that ioflsence ®f these strains on the tempe- 
rature itself it not tee large. Thss osississ of sncbv** 
nicol coupling tares is energy equation weald be messing 
less* 


second simplifying essumption is omission of 
inertia term in m equation of motioo. Duheael had 
shew** that the time rote of change of temperatnre |s 
slow enough so that this term should sot bo signifies* t* 
Danilevskaya ^ end Beley have also shewed that ©mi- 
ssion of iwortin term is meaningless The two simp 11* 
fyiog as sumptions give rise to nncenpled quasi-static 
theory* This on coupled qeesi-stetie theory is divided 
lots' two di stine t subjects, sad known os theories of heat 
conducting end the rmoelasti city. So lotion of host eon* 
due ties aquation gives the distribution of temper stare 
I welds the body. 

la the present esse, it is assameti that leaps* 
rota re distribution in side the thin shallow shall is 
known before bead* It is tsken that temper a taro is ssti* 
symmetric cheat the middle plena of the thta shallow 
shell bet otherwise arbitrary* 

itoa to elastic cess teats each on Young's SSedelms, 


? 

fho tbemal stresses era related in a cosap Heeled 
nonaer to the temperature distribution. Blech an ic el 
end Metelltirgieol properties ®sy bo changed with the 
tea per store. Xt is assumed that eleetle cone tea to 
do not depend upon the temperature m& notorial behave s 
elastically at ell Uat* 

only Gradewcsy**^ hae tried to derive the 
governing equations for thin shallow shells. Ha hoe 
tehee Kaviosr type solution for da fl action of thin shallow 
shell* Kevlar type solution Is possible only for all 
edges simply supported. Set Mavlor type soletloe is is 
terras of doable trigeo oac trie series which believe* poor* 
ly near the end poist of interval* Cradetrczyfe ha* set 
salved eosplotly for ay type of boundary conditions* 

Stress function end deflection of thin shallow 
shell subjected to t tape rater* seat satisfy the following 
two equations. 

&1 V s V s w- t* c 9 ) « ^ - V s J ...*<8) 

and 

* g |» L 0») * V s Wf . .(3) 

Stress and deflection functions are separated 
into different equations. Solution* of these equations 
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are obtained ns t.avy’s type solutions which ere in 
t«nu of hyperbolic s»d trlgono©© trie fanctio-na. It 
gives hotter convergence sear the end la comparison 
to Wavier type solutions which is in l arias of double 
tr igneous trie series. 

Ibe boundary conditions considered or#} 

a. two opposite edge* of thia shell ova shell simply 
supported md remaining two edges supper ted io 
my &m eer, 

b. All four edges of thia shallow shell simply 
supported, 

e« Hie opposite edges of this shallow shell simply 

supported end remaining two edges c looped. 
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a i/4>m n 

gamuLfeiiwi or besic ecu axioms 

Following assertions Nave been made la the 
derivation of basic equations; 

1* Tbs material is homegenosast isotropic aad 

obeys Heok*s LiMi 

f>. Planes peraeadl color to middle sarfaee before 

bending* remain plane and perpendicular to middle sur- 
face after bending end de net stretch in their length* 
l.e. <*3*0 , i S3 « 0, G s * 0 

3. WOnael stresses on planes parallel to middle 

surface ere small in coaporisien to other stresses 
and these normal stresses ere neglected. 

<*3 m 0 

4* inflections are smell compared te thickness sad 

the effect of loads in the piano of shell is neglected 
and coordinates of a particle before sod otter de- 
fa motion are the same. 

6. The elastic constants Strand the coefficient 

of t hemal expansion ore independent of tempo retard. 

6. The thieknoss of the thin shell fte very smell 

in compsrision with the least radios of eorvetore of 
the middle surface. 

i.e. 1 * V*i • 1 and 1 1 V R & * 1 
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In fchwi above 3;s§> region, equating the eoeffi ©lento of 
6 Uq* 6«q ®md 6w equal to koxo, gives bho equilibrium 
equations* 



bA s 

0 ?, 


*2 + 



S + 


0 

**•*■ ( A|S) + A^A^q^ V 0 
S» 

IlMHIl ( 16 ) 


ft 8 A^ 8 A 5 & 

6v o * (^a J * ,7 *1*76, Ca 2 s > * A i*a «g ” 9 


1 » 1» 1 


»«**»«*•« 


( 17 ) 


&wt Si| 


^v s «+!l+! 2 - w-v 2 !* 

8^*2 ( 1 «u> ) 


**••«•«• * 


(* 3 ) 
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where. 


v 2 (...) 



vanishing of surface integral term© in elation (1*5) give the 
appropriate boundry conditions considering the nurfe.ee integrals, 
to get the conditions ©Ion:? ^ « 0 end ^ , the surface 

integral cmat vanish. 

/ ^ 2 {^ t 6tt 0 ♦ 3 :'v 0 - jJ ||^ * (U f f->*}Wg » 0 (l<Wa) 


2 This can be done by obcosing 


El $n e «■ 0 } S r >v 0 * 0 (19-to) 


where 



As celled Kirahheff'© abser* 


This physically m®m that the virtual m$k don® by farce® 
end moment® ©long the bouadrie© » 0 end pj * gj J®uat be aero# 
Thin ecu be done in three ways, firstly by aaklsg all forces end 
moment© to be aero along the boundry or aoooatfly by choosing the 
diopleoemeate and slope© along the boundry to be sere* or 'thirdly 
by choosing partly forces and moms at© to be mm and portly dis- 
placements and slopes to be aero along the boundry. 


Similarly vanishing of surface integral for m 0 and 
P a *' 1© given 
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/ A,f 5j r»T 0 + s -•& |jj» 


1 ®»49 

(B g + i*- — ^ ) 5* } d;^ » 0 


(20-®) 


Boundary conditions alon# ?,g • 0 and pg ** con b#; obtained aa 




(20-b) 


At the few comers# 

2K<|g ® 0 


(20»e) 


Consider a ©tree® function *9 T which l» rotated to the 
lnplene forces m follow® i 



tutting the ueluea of T g and * into «tp©* ( 16 )# ( 17 ) 
end ( 16 ) end simplifying by uaiuf the ^tjee-Codassi relation®^ 



(22) 

( 25 ) 
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In eboenoe of forces in end f; g directions* £qn** (S3) 
and (23) will b© automatically satisfied for shell of mm 
Gnu® slan curvature* for shell of non »*ro gauaalea curvature* 
the ©hove equations are approximately ©atiefiod for shallow shell 0 * 
In order to obtain* the unknown stress function , o , # am? w® have 
to consider the ©osapatibility conditions. 

Compatibility Equation to given by 


♦ T*q C 4; ^ { *8 wf * ^ (ej - e°> . 

, ee®, a*, « * ee® 

~r A i e tsl* rjfj -q i 4 nrjj- * 

1f~ < e t - *a> - I *t "iff - wf %]] * 0 (0S) 

By substituting the value© of fc g# * S g * €® g in Eqn. (25)# 

"*<■» ♦ v8 **» ♦ JJSf L > ♦ 

- A ■ 

|f 1) . J L IlLi^l^ 1i 1 M t 5*15 ^ 

TjJ { Tp% 1J *tj }J % “■'. 

for * 0*9 Gmuoalaa curvature' (l/HjBg » 0) 

. «&(»)' ♦ v a v^ » v 2 % (tr) 


(These equations hoi cl good tor 

1. shells of «-ro G. r uesian ourvtf-uzs swob as cylindrical 
shells and cones, 

2, good for shallow shells 

5, approximate for non-sera "mmvim curvature* 

These ©nuatioao aw coupled carl requires the eicultencoua 
solution of deflection *w* and stress fu ction If there 

is no external loading*, for constant thickness* equation (24) 
reduces to the fores* 

J>,| -» L ( 0 ) ** •» (!?8) 

Superposition of deflections due to transverse lead alone and 
deflection du® to temperature distribution alone && possible. 
Solution of deflection duo to tenperstuxe d is trihut ion ©Ion® 
has been utilised in the present case. 


CHAFER III 


tfifliiOB Off JOluaioas 

In thin chapter* solutions to obtain thorrol strokes 
in thin shallow shells with two •© d&an ainpl.v artf other two 
edges supported in ssy manner are d«wio?Hii, Th® equation® 
gov raing, the t ran sverme dlfleotioa of elastic thin ©hello 
subjected to ® temperature distribution which is anti-spametric 


©bout the middle plane and otherwise arbitrary ares 

B,| V 2 V*w *• !(©) ♦ m 0 (29) 

lh£(w) + V 2 V 2 © - V 2 % * 0 (SO) 

Be fleet ion function *ir* and et retro fUnotlcn *0* may be 
mp ® rated out ms 

D t v e © * m 2 (e) * % ♦ mv 2 1 ) (SI) 

© 1 V 0 W * Etol 2 <w) • ~^7 6 ■flJvT ♦ vh €%) (32) 


Both functions w mad © « hew the saw Imnogemsteus ©olutiom 
but differing in constants of integration©, 

Doubly curved thin shells may to* raduo&l to thin Shallow 
shells* when It is satisfying the following conditions 

$ * s(x t y) « M < «J* (m* to) 

<&* «> , t#> 2 « 1 , < H - # > « i 


Thin ©hallow shell hem b*ea ©he«n la fi-?* (8) with 


coord larto »y©t«n* !i'hs value of rnd for ©hwllaw shall 
gas>,y ho approximated as 

/if «w Ag ® t# mil Rg *rs eonatent* 

iitJttae couples and atrtas rcsultrat© for thin ©hallow 
shall in tortaa of do fleet Ion •w* find otrsoa function Q ssay b© 
exoroeASd o® 


iL, » 

** 

(£fs 

e * 2 

2 ^, 

03) 

1 * m 

y 

+ 33 1 

<i!a 4. v 

w 

# 2 ?K % 

Os) 

M „ nw 

ny 

* a,(i 



( 35 ) 

T 1 * 

i!a 

•» *>Hg 

• *8 ” 

ft®. A ., g®fl 

•y. Jj $& ■WWflSw 

( 56 -n-e) 

ey" 

&x 4 ’ 6x^y 


Th? middle plan® ©train® is terns of 
he expressed ns, 

®1 - A < *1 " V *2 + V » £ ? " 


(Ifg *^9? | * %) 


Xtoubly curved thin shallow 
s « 0 end x * a and waking ^ 


:.(?)- o «aa eg (•) - o 





or 

e<£. y) - o (30) 

distribution of teaoeriv to.ro iw ^tii^Erairie ebou* tb© 


micHI 

l.c plea© eo that ‘1. 

O - 0 



*.o /o „ e c o 

W, ( r »y) * — 4y 

i‘-!X 

{ * 0 
;'JC «• 0,0 


or 

_£« (I . y) . 

fcas* 

0 

(30) 


w •* Ot Mjj « 0 

when x * o* & 

(40) 


Boundary oonditious (40) imply that 

» » oj a-5 m « JMo*y)/&i(l~ y ) fttoaft a* o 
ear 1 

S ^|U 

®»0| iiioa^ 3? » m 

2he boundary condition© ore aoiwhoBsogcnooua. The ac*>-ho«og*no ity 
is xcnovftd by tskiag th© Amotion *»* into two itonas^®* st 
follows 

»(** y) » u(*» y)-+ H^x) V 2^(e,y) ft ft (aO^ (41) 

0 -v> O-v) 


Applying tbs boundary conditions into ths equation (41) 


E(0* y) « * 
y) 


[MjCO, y) H a (0) * t^(a , y> H,(0)] (48-e) 

0-v) C'-v> 

w - % (0, y) V®> 

oTvo cT-v; 

(42J&) 
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ll(0 # y) * ~[^(o,y) H 0 (&) * <^{Q»y) J* ft (n)] (42-n> 

Cr-V) O-v) 

’"“’“R ^ ® ** ^4r (0»y) ^15 *• Rr £ &f$0 ** U«( © 

e* 4 - <br 0 f— "1 a®*' 

C/-VJ 

» * * {4<wi} 

Function® R 0 (x) end H ft (x) are chosen in the following iaf.ar*eaf Its 
order to nokr the equation (42) horwjteneoua* 


H o (0) » D 

h b (0) - 0 

*4j H o (0) * - 1/^ 

dsc 

l '» (0> 

(43-e) 


KL(a) » 0 

Q 

'^(a) » 0 





H q (©) ** 0 


2 

-4f « 1/6« 

dx* 


She fun ot long H (x) and H fe (x) are taken ae nolynonicle 
mtth four constants which are dot am" it® 3 tilth Blow condition® 


on U & (x) and B e (ac), *h«y are of the tom 



(j 44**^0 

(44~b) 


the boundary condition® on t?(s, y) b«ocn», 

U « o t « 0 along ae » 0 and sc «* 0, 

' lx® 



Now substituting the eq.n. (41) to (32) 


where 


» 1 v 8 ^ ♦ Ih& 2 («) • F <x,y) 

F(x,y) * -V 6 + V a L(^) 

(t«v) 


(45~*0 

r i^ 6 C I! c|[Ot»)8 0 (*> 

0=0 


* %(«»y) H e (x)] ~ '^(x) 

(pv; 77TP; 

+ ^(&»y) H e (x)j {4S~h> 

o 1 ^) 

Solution of ©qn* (45**©) Key br teV.oa ?•© 

U(x,y) ® z T ra (y) sin ~jp (46) 


Ibis form of ©olutlcn will m Ai»£y the bountaxy ©osutl- 
ti c,m ©long x ® 0 and x » a* 

Substituting tte wive of 0(x| y) in the equation 43(») 
©Oil Btoplifyisi^ 


t ■$ V y) - ^ v»> > 4 

Ih a 2 - -JS* 2 « 

"7 «£:!* V*>~ x V*> > 8 *V*> <**•> 



(47^>) 


/:® 0 uss@ ^(y) in the form 


Y (y) m o 
M vjr/ K 


tobetituting velue of Y m (y) into the elation (47«*e).$fvs& 


1m 

,B 

WB 

+ ^2n + 

•;* * 4$k * 

a, 3« ,? ifi * 

3 te 

sn. & 

w » 

* b 5b “ 1 

l 1m 

m 

1 


A** 

*a 

«r 

B 1a 

2ia 

m 

-4 AA 


A 2m 

ia? 

V 8 

‘3® 

m 

(6 m 4 n 4 + T 

) 

h 3® 

a» 

B V ft4 


m 

~2AA(2bAi 

4 ♦ ?E) 

A 4b 

m 

B.^/a 6 

*5* 

m 


■ R%*) 


m 

W a ® 


f * 12( 1-V ) 2 . H «* K 4 /IU 


thi® equation into fourth order equation and 
method given in Ref. (9)« 


solving by the 


Depending vmikyalue® of A^ m to eqn* (47-o) hea 
different set of roots and different type of function Y R (y) as 
given in the Appendix (A), 

the homogeneous eolution^ 1 ^ Of ®qn. (49) ney be espreoood 




alon& y ® o 


rVjp it; 

The superset 


[oj » 


Tfe© differ®: 

i. 

[%] 
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lor hoROgeneous howndery conditions ©t at * 0, end x ** a 


0 (x»y) ® 


Off* 


j 3 


? m *(y) ein utat/a 


T 


2P R (y) * [a] [* ri (y)“ 


(54) 

(55) 


[B] 





Thera ate sixteen unknown®* eight in deflection *w* end eight 
in stress function *e». Bat possible nucber of ho nndery 
conditions art eight along the edges |» 0 ana y ®» b„ Eight tmlmoim® 
©an he expressed in ttra ©f ether eight nateacwns m giwn in 
reference <tt), 

Substitute the ’eel no of 8(x*y) ®s»& hosiegemoue boundary 
condition pert of w( 3 c,y) in the egmticn (30) 

^ ~"g?~ * ^ C*bJ 

>f^# *i%i 4 - ^ihi 2 }m f L^] m ^ 

-..(so 
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;,'o ©atlefy this © quation, eguel coefficient of y^ r 
through y^ equel to aero, ’ibis giynn eight equation® relating 
ttee eoaa:eate c t B * c 2n* J 3n» C 4»* S&ra* a 7T?* ’"'era to 

%m> J: 4m f *W %»• V %m- r W*wta« ttj * B 1© % 

coaotcnt* into 0^ o through 0^, T-=ua through arc 3motm 

inters of 04,3 through c^. fad there ?w eirfht uafenenm 
through Cg :n . 

g.quation (47-b) cs a bo written r.s 


(iMkj (B-k 2 ) (Mg) *» E{y) (57) 

whore a^, ® s or® the ro&ia of the aurlelly equation 

Put u 1 ■ (D-&g) ... (XU5 a ) T a (y) 

Bo that equation (57) become® 

(D-*^) U 1 w K(y) 

This ir. e fin border linear equation in the unknown for 
which a particular integral ie xeadly found, i?ejrt put 

Ug « (IM^) (B*a 4 ) ,* (5J-ftg) \j(y) 

and write 

(JV»g) Ug * 

where i© known funotion of y, ami mlm this ft rot order 
equation for Ug continuing la tbie a. ©equenm of linear first 
order equation© are obtained* 

(i * i t g| ...a) (53) 





vhex® 

u c * £(y) 8 Up » 

....(I^c e ) ^(y) (A =* 1,S, ...7) 

and v-her©, for ©cob wlue of 1, 1 r know*? function of y. 

/..ft 't 1 eewn each etepo, obtria the fir* ’ f* ret e^lcr differential 
equt-ticn of the sequence (58) » ne»eX.v 

(r-e. 8 ) \Cy> - « 7 

a solution of which is a particular ir!tr-.;-;rol of (5Y}» The 
sequential integration of the ©qn&tisn (58) ia obviously a 


routine process. I'aoh ouch equation fcre the integrating factor 


a 

i 


i* 


and the solution 


» e 


~&*y 

/ • ” u^Cy) dy (i « 1* ...S) 


The unknown constant® through any bo obtained by applying 
the suitable eight initial perimeters* diopleoetsents, slopes, force®, 
and ssciBsents in tews of Yg,<y), * \ %1 (y)* ^ a (y)*3P 0 1 <S r ) 

and F 0 ^(y) along the edges y * 0 and y » 6, depending upon the 
boundary condition® ot y » 0 and y « S, 


This may bo expressed in B&trix foxsn 


[OJ 


■ U] Lcj 
» [A ] -1 


[*!»] 


(59) 

(« 0 ) 


The t®ln© of t nm mtTix dsptMsj upon fhe boundary 

boon gimst for simply 
y » 0 and 


condition®, V&luos of 1.^ aateix 


supported along y « 0 and y » b 

y « 6* 



Sfctrix CA] is given as 


Vi 

6 

IP 

A,| 

14 

A 15 


V? 


V i 

A i>2 

A 23 

Ar> , 

r-** 

Vs 

Vs 

V? 

tL& 

Vi 

A** ®*> 


ft 

-'34 

*35 

M , 

•A r rtf 1 

•37 

* 

Vi 

I 42 

A 43 

A a 

Vs 

^ ' jl •*" 

k !y> 

V? 

Va 

Vi 

Va 

Vs 

V 4 

Vs 

Ay. 

' 7" / 

V? 

Va 

Vi 

%2 

Va 

V 4 

Vs 

ft 

Vv 

%f. 

Vi 

Vp 

Vs 

A 74 

Vs 

if J .'4'!' 

iy 

Vt 

A 7B 

Vi 

V? 

Vs 

V* 

V,*** 

f t-P 

V? 

Ve 


Velue of Matrix til denendw ucv,n[l£p(p )] %pe of boundary 
conditions ©long y • 0 and y * ®. 

Value© of eleseutfi of Matrix [d-J have boon, given to 
txanple for ®i«ply arup ported alone y®0 and y * b. 

Hut the sake of aiwylioity* *sa&- teaporetutv distribu- 
tion independent of x and y Eqn, (4?-«b} reduce® to borogemous 
fora. S© deflection function aay b© given by 

®» 1 $. 

w(x, y) «■ 2 (" T m (y) * 4 —■» — A *" ■■■ ] tie ®sc/« 

»* 1 B,(1-v)»V 



the boundary condition® ere* 
w(x» £ ) * 0« 


fbU fives T te <©) 


t«,Cb) » 
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Botmd&Ty aoaAfctioJMi !*re 


*{*» §) «-0 
•IM® gives t ffi (0) 


4 VL »® 

<* JSftp ® @'Sd s&Slly 



® 0 


Shle fsiwa ^ ■yo) « ^ Y n <b) 

* a (*« £) ■ © » ■thin Kiffe f R ( jj ) » 0 

U(X* J) * 0 , tMs g!v« -V,(S) m 0. 

ay 

- 4*?^ ir 

- 4 ^ •^/EjCl-'OiAt* 


f V 3 ~ 


o 

0 

0 

0 

0 

0 


T 


* C °3 LV y) ] 


m 

% 

a»1 


4 1^, *? 


( "S> *** «*\ 

Y«(y> + J&tsc/a (61 > 

9 B,(W1 ®rV J 


" fjD *-§- Co] T + Da ] 1 [K R ] 2 }(r c (y')] Bias 
'*•> 18 ' “ (6C) 

x ( #, {- lo) ' ty 2 ♦ ^ ■[«] * } [?«<?>] 



4t-v)IW 


aia 


{655 
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t. [ 3 j £ .»Wj CO 3 BlT!X/a (64) 

( a) ffseKnle s 

fry 

In order to illustrate the n*t.h©-:! # en exnsrjlc v.itb © fc '/4H^*5 
v *# ,5, end I^/8g • 0*5 hm 'been v/ori'e-i out fasten • 

Hoots of the equation (49-a) featv the io.rr. «ig ■«• ♦ n g 

t (n^ + in^) 

Y ?i (y) * ®£cosh n^y/e, ©in n^y/e, nosh n^y/e, ©inh a^jr/e, 
cosfe n^y/n cos n^y/w $ a inti n^y/a 00 © n^y/e, 

©osb n^y/u sin a^y/e, ©into n^y/e ©in n^y/a } 
Fqimtlng the ooeff iciest© of elerenta of t r>tris [ Y«(y)]^ in 


©fillet ion (56) to sero» 

constant© 

D,|. r , to ©re obtained interms 

°* 8 1c t0 :: 8m " 




c 1« " ® 9 

^1* 

* ®2o ’ 

®3 °2* 

“ ®10 

c 5® 

• B 4» “ 

K, 'l0 S 4B 


®4m m K < p ?%> + P S W • %« w k < p 7 <fc* * % V 

*?» 49 «*3 %» * % V» *W - *<**Si ♦ > 

P| » - Rf/fig a?« S + n| - »| 

*2 - *»,» 4 

t 

Fj « aV ♦ a| ♦ »J •* Snjj - ®V <n| ** si| ) 

*4 « 4n| a 4 - 4»J aj - 2 A® 


where 
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i',- ** ^ / * 

i’g as ( + P^P^)/ 1 * 

P ? - + p 2 p 4 >/ p 

% - <V* * P 1 P 4 >/ P . 

Vq ® (•H.j/Rg m2 ^ 8 * tfj) / + qJ * R®tt® » P ) 

^10 * ® 2rS ♦ n|) / (®*n* *■ n| - ttfr? ni) 


&h 


K « 


p £2 

? - Pf v *4 


totals [A] ha© different valuer for different typaef of boundary 


conditions. 


®lori ' y ® 0 fsnd y » b 


■^IK * ^ < 


~13 ” -15 
Coah n^b/a 


*25 “ 


©ash a^b/a 

eo®h n^b/a coi n^b/a 
000 Rgb/a ©in a^b/a 


l 18 " *14 " ■*« “ *17 “ *1 


Agp » Slab 


$2$ ® Slab a. 


JW 43 


®&WB co® a, 


‘\n <® slab a^b/a oia a^b/a 


^31 a 


m 2 „2 

°3 * a 4 


A 3S “ 


Hi a 


Hs * 


a i A ai 

a i *23 


^ ^ «o 

i^#*’ & 


« aj #23* 


*2 *S>4 
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*46 

ra 

A 35 

% * p £ 

1 A ?7 

if> 

l m 

S3 

'?5 

*ea - 

1*1 

•"'2 

*2B 

'47 

S3 

H* 

t P ? 

h& 

V' 

33 

; 35 


•fj 

2 

1 "25 

Hi 

m 

h 

* 



%3 

es 

"Ho 




A 5S 

BS 

p 7 




*50 


P 8 




A 58 

f» 

Ha 

» 0 



A 56 

«a 

A §7 

©J $ 




ISt 

"9 

Vi 




®Sf 

V 

*22 




$8 

p 10 

V? 



Ha 

83 

*10 




%5 

m 

■p 

’7 

hs * 

P 5 

Vs 


*» 

V 

*26 * 

*5 

hi 

Hi 

18 

P 8 

A 26 * 

p 6 

hi 

Hn 

53 

V 

*25 * 

% 

A £8 

M n 

cs 

*9 




% 

33 

P 10 

2 

a i 



a 75 

tS3 

x> 7 

k 35 * 


s 

*78 

ts 

% ’ 

'•55 ! 

P«T 

<$£ 

6 

A 72 

m 

A 74 

w 0 



% 

«3# 


t* fj 



i- ©1 

8SJ 

*9 

A 41 



*32 


F 9 




% 

Ml 

P 1Q 

i A 43 



*84 

s» 

P 10 

*44 



% 

Ml 

P 7 

A 45 * 

P 5 

A 46 

% 

» 

P ? 

% * 

? 9 

*46 


63 

P 8 

A 46 * 

P 6 

Hi 

*$S 

«» 

% 

*45 "* 

p * 

* 4 ® 


Solution of the ©atria equities* <60), by pitting th© values 
of matrix C *p m l matrix [a 3 , glw» the vb%wb of constants* 

C 1ffl through c^ # flsilootioa ©ad ©©swat® say fe© obtained by sub®- 
tltutlcm of 0^ through <3^ .in Bep« <#?).# <68 t <6f) and (&&)« 
Value© of olestshh of the matrix Cdl* for shallow ©hell with two 
edge® ©imply supported end rosaintog l» © d#$©t ©lathed, hav® boon 
.gitwu in oosputer pregeai®® in Appendix 0, 


r ihe ©etbod Jc 


is toot# rrnerrO. then /.rbartau ry&H’f* 


method* It does not re iuire to define deflection la tea® of 
©ny other functions such eg w » V\) C) Where is any 

arbitrary function* Stress function ®e* end dr flection function 
•w* err replaced in terras of ty ftenotitcj In F^* (£#) in rts&ert** 
a tray an* a method, 


*5 toe- solution of differential equation which is in terra 
of function, i*» obtained with the help of boundary coadltieao. 
Then stress function sad deflection *v/ 4 are obtain® d froia 
function. But the solution of differential equation which is in 
teres of 4 > function is test henat loclly diffiault or reont pro** 
bcbly iaposrlble for cones Where SL is not sero* Method 
developed in proceeding actions isny be applied for cocos whero 


method hold® good except for addition of particular integral 
term in Iq, (50)* 

p . 

Parameter & /h ft* end b/a ratio play -cm important role 
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two edges b imply supported e-ad resptnin;? t>. o edfjos clscof-d • 
ifrilerly* for fixed value of parameter r^/h diasensicnless 
quantity for di fleet ion of the shallow shell fro© 

upward direction to down ward direction ©s values of the b/& 
ratio increases as shown in. Appendix r * for oil edges simply 
supported end two edges steely supported end renpinine; iwooigoe 
clamped. Values of dimexmionloss quantities for end h 
decreases in Magnitude with increasing the parameter a 2 /h P, 
end b/a ratio* But change in Magnitude due to notesseter 
©ad b/a ratio of dimensionless® quantities of nom&ts B x and IL 
is ©mil in comparison to the ohenge in magnitude rail direction 
of dimensionless quantity for deflect ion* Magnitude end direc- 


tion of the dimensionless quantity for moment s?^ also change with 
change of parameter * 2 /h s. »«« V ratio as nhown in Appendix ®* 
For all caeca numerical %lue© and graph® have hvm pMseatsd for 


deflection and Bonents (M % and Ky) at x * e/2 w! 1^/Rg * 0.9 for vari 
various mime® of parameter a®/hB 4 and b/n ratio* 5Taie has been 
done for parameter a 2 /hF 1 ranging feasi 9*0 to ro,0 *t inteJml® 
of 5*0* and t»/a ratio arranging from 1*0 to 2,0 at interval of 0*2* 

S’h© deflection end nossnts hav® been presented by the diwaoion- 
less peraneters wb/affa 2 « tL/feafh 2 * U-JlZtfih? m& JBU/Baffc 2 , 
Ilmnsloalsan >•»*»» V Mb* baa U« ««« * » 

Only magnitude of dlmxmienles* parameter® t^/tecSfe » tfy/Fidbr 
have been presented in Appendix E and direction in given in • 
arable Appendix S« Since the Variation of ub/effs** J^/BaSto 2 * 

My/f ofh^ along x-@sle is in simple etxn> form, the values of then® 
Quantities at any point can bo easily circulated. 



She r©0al,tfi obtained in prooeodlrt."? e eoti^-Jin, caay toe used 
la desire of shallow shell* e roof a whore* outside the building* 
etsos ;herlo terneretture le different fyom inside tetsnereiiar 
of the building. She aethod developed In preceedixi# eectione 
Key tor weed .in finding out the itreonet in end cepe of tooilere 
end eirwing® where temperature fircdicat estate, duo to tem;jerr.«. 
turf; difference of inside and outside surfaces 
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fu n^miK a 


CHA^/.CiTliKIwTXC KOG?,: 0? TKi; ■ fr'-SSE:/: IAL '^.QA'-T:.:; (47) 


i>1 ite fcoKOgeneoua part of the Atffcrcntiel ©•••’ur.tion (47-h) oca 

bo ^jittPW vn 


»0 J> .A $** 


V»> " ; -3K V y) ' 


* a 5i: \<y) * o 


U-1> 


where;- 

* to ar ® $iv©n by eqa, (4%b) 

The nature of the eight oherrotortetio root* of rq# (iut) 
depend o on the value of to ©n4 can be daeslfled under 
th® felloe ini; ease*. 


Cm©© 1 

+ »r 

i*V 

+ 

a 3* 

± 

n 4 

C&a© 2 

+ 

* n 2» 

* 

*V 


la 

Owe 3 

± 

1 V 

* 

hij» 

Hh 

In 

i 


Case 4 ♦ a^# ± t»2» * in^* ± i»| 

Ceee 5 ♦ 1®I# + iUg t ± tn^» * Aa^ 

Case 6 + n^» + »g» ♦ {n^ ± in^) 

Gee© 7 ± n^ 9 + lttg, ± (n^ ± in 4 ) 

Owe © + ia^t ♦ i»g 9 ± (n^ * in^) 

Cage s + (n^ * fa^)# ± * ta^> 

whose n^, Hg # end a^ ere the oo«i>o&i*at* of ttm oharaotariotlo 


root®* 



t.2 , 1 jilt* It *. 3? 
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?h« baste of Rolut Ion® [¥yj(y!iJ * (*>?. (91-e) } t« given 
be lav as a row ©atria E*r,(y)] * , eloi.v--.itfc the set*!* [-^ _ 
(Kq* (52)5 > for Ti-riouB eaten Hate* 1 in *’ene, (-"«") 


Ciaso 1 



^coefc n^y 
oosh n^y 


slab , ccrh n^y, sinh n^y 

einfc n«y , co^h n^y, sinh (/~3«) 


o 



0 

0 

0 

0 

0 

0 


d f n 1 OOOOOO 

0 0 0 0^00 
0 0 fi 0 JU 0 0 0 0 

0 Up 0 0 0 C? 0 

0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 f^a 4 

0 0 0 0 0 a* 0 

4 -l 


<*~3b) 


where 




[■5f H (y)J * * {oo»h » t y, 
coefc n^y» 


slab n^, coeh f einh a^y, 

sinh n^y* oo® n^y , sin n^yi {Mft) 


QWjJ is given by Eq, (4-3b) with 


4^ » dg » 4j "» 1 sad * «*1 


( A-4h) 
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'.;®se 3 

[Vy>J * 


jcosh n^y , Bin':- a^y , Booh svj 
cow a^y * sin a^y * oo« n^y 


cinh n^y* 
sin n^y 


( ) 

given by Eq. (/-5b) wiV? 

* ^2 « 1 » frtttf 4 ^ ** <tj *»1 (/-.^b) 

Case •; 

[^.(y/] r • ^oo*h » t y » slab x*^ , can n ? y » aim a^y , 

co® n^y » ©in a^y , cos a^y s ala n ^} ( iU4fe) 

is glvf-a byi £q* (Mb) with 

« In aai dg « « -1 (^«6by 


Case 5 

[■y R (y)] S * on nfi 4 «ln n^y t »?» a-,y 
o;« a-»y * ala n^y , coa a^y 

£ ie given by ttu (JU3b) with 

*1*" a a * *5 * a 4 E ' - 1 


ain a.-y 

4 - 

ala 7. (JU7&) 




©roan 


x 

[5 R (y)] ® ^cosb n^y , «iah a^y .» , ^ 

ooafe a^y oo* a^y > eiah co® a^y 

©lab n^y . ©in a 4 y * cash ji^y aln a 4 y](A~€ 


Case 6 
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0 

n \ 

0 




0 

0 

0 


o 


0 0 0 0 



0 


0 


0 




u 


0 a,, 

0 o 

0 o 


0 0 

f\ S'l 

• J 


0 


*5 


0 0 0 0 

0 0 0 9, 

V 


■3 0 

0 o 


C 0 

0 * r ~ ! l ,n 

■ > 





•‘2 


C ^«8b) 


whers 

a 1 “ *8 " 1 

C&©o 7 

[^(y)] ^ » {cosh a^y , ©into a^y v oos »<-■ , win n^y f 

cosh a^y > eos n^y * ©Sab n^y, cor a^y » 

»inh a-jr , sin si^y » oouh n^y, eln a^y ]. (*-%) 



fi ® 1, ©ail dg « -1 (vft#9b> 


Oft 86 B 

[t K (y>] * * £«oe a^y # sin n^y , coa n g y , ftln n^y* 

oosb a^y oos? n,y „ sin »^y cost n^y 9 

®lnh n^y oln n^y # cosh n^y Bin »^yj (A«lOa) 


A* fives by Iq» (Mb) with 

&1 « dg * «»i . 
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C®»# 9 

£* R (jr)] T *|* cosh ttjjr sos agr, slab ujy 
eosh 8jf sia Bgjp, slab a^y 
cosh a^y cos k^ t stab a^y 
«Mh a^y ala a^ ( slab ay 



COS krf, 

sia 

C®8 djja 

aia a^y } 



(A-lla) 


(A- life) 



»?mxx % 


FLOs* CHARI'S 


ftfoia Flow Chart 
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§_§?£! 2 kY_su i j . r n a c s c l r c r. l i s r 

' r;T: T ” T ~~ PM 

I R r T C COLFCK 

r.' i r-* v f . £ .; i rj;i c, r r w . p t ( a ) , Ar-t a j .p.oia) ,x i (A j , ac ( '♦) ,xr(4 ) 

:: I J'Ff fiTCU TUFT A ) » HP (4 > ,KR (4 ) ,Fk( 4, 2 ) t F l (4 ,2) ,PP( A } 

:. i x f !. ’■ [ r,:> I r* j vcT (?) , a l n , ft ) , u 1 3 , 1 ) , I *i Of-x I0,2> ,Plvnn*> 

r I >'!. r sir ;i -M 1 0 ) , ’,< ( 9 ) » P { S ) , Xf' { ‘5 ) , 0( 4 ) ♦ <vr-H 7 ) 

F- ;/ .'i 

= i . c 

? J K Ii-M < 

3 rr 7 f i* j 1 > 

4 i> - 0. 0 

;< - o . c 
? p = l . n 

7 « l> I! = 3 . !. 4 1 6 
FACT!;=t»I I* PM 

v=n. 7 

'.= 12.-- { l.Q-V**Z)*X**2 

0 ( 1 I = 1 . c 

C(2)=-4.I)*U-ACTR ) **2 
C [ 3 ) - 1 . 0 * ( F A C.TR ) * * 4 + L 

C(4 )=-?♦* (FACTk )**?*( (2. *FACTR )**4+U*R ) 

C ( 5 ) = { F AC TR )**4*((FACTR )**4+R**2*U) 

A 3 =C ( 2 ) /C ( 1 ) 

A2= r .(3)/C(U 

ai-guj/cu) 

IKIil/CU) 

AT = A 3 / 2 . 

Am )-i. 

AC f ?)=-A2 
AC( 3> = A1*A3-4.*A(7 
AC ( 4 ) =A0* ( 4 , * A 2- A 3**2 }-A 1 **2 
CALL CUL'IC (AC,KT«RT1) 

IF (RTF >20,10, 20 
LC IF{«Hl)-RT(2))ll t 12,12 
1) RT { 1 ) =RT C 2 ) 

12 IF (R T ( 1 J— k T C 3 I >13,20,20 

n RTm=rtT{3> 

20 C=RT(l)/2. 

IF (0**2-401 22, 22, 24 
22 C-C* 

CA-SGRT ( AT**2 + 2.*C-A? > 
i-C Tfj 25 

24 C=SCRT (0**2-ACJ 
C A =- ( A l / ? .— A T * C ) / n 

2c 4cm = i. 

AQ ( 2 ) =AT-C A 
AG ( 3 ) =0-D 

CALL 0/40(40, XRI1 >,XR(2) *XI M.) > 

DC m=i . 

RC ( 2 ) =AT+C A , 

RC(3)=G+D 

CALL CUAD(BO»XR( 3 > ,XR(4),XI(3>) 

XH2)=-XI 1 1 } 

XI 14)^X113) 

PRINT 7 7 , X : 

77 FORMAT ( IX ,,F7 . 2 ) 
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te- 


T ' r !'f" :-MT 


F H K T r < A ?J s i'IJRC E LIST 


l- S 

FOP. ‘-'A T (1 X ,F7. 2 } 

COL 1=1,4 

pim I) = < l ( 1 )/;<>. [ i ) 

T h F T A C I )»AT«MFP( I ) ) 

PI =/,.*.*, TAN( 1 . ) 

( l ) =S'j’-! T ( X? { ! )**2 + xl < I ) **2 ) 

O? { I ] =Si,P T ( OP ( I ) ) 

f P ( ( , I ) = C I t *Cf j5 (Pl+0. ‘i* THE T A (HI 
Fin ,21=8.01 I )*SIM2.*PI+0.5*THETA{ I ) ) 

FOC ! I )*Cnr>t2.*ni + n.5*THETA( I ) 1 

J OTP , I l F;l U , J I , J = 1 , t , ) , J = 1 , 2 ) 

!>:; iwt a 1 1 ( r- 1 u , j ) , i * i , 't ) , j= i , z j 

Ff.'iWATt IX .PE 14. 4) 

IMTIALIZATICf, 

(A IJsHU ,?> 

'A?) = F?l ( 2 ,?) 
r :C3) = Fk n,2) 
ntA )=r-i f/ + ,? ) 

0=3. 1416 

P(l)=-P *(PM#G)**2+ Q(3)«2- Q(4)**2 
Pt2»=2-*0(3)*C(4) 

P(3) = (PP*G)**4+Q[3)**4 + '3l4)**4-6.*( 0( 31*0(4) ) **2- 1 Q ( 3 )**2-L) ( 4 } **2 ) 

2 * (PM* 0**2 

P ( A ) = 4 . *0 l 3 )**3*UI 4)-4 .*G(4)**3*Ql .3 )-2.*Q| 4)*(j( 3)*{ PK*G)**2 

PK = P f 3 ) **2+P < 4 ) #*2 

P(5) = (Pm*P(ll-P*3)*Pl2) l/PK 

PC A 1 = (IMP )*p (4 1 +P< U*F ( 3)1 /PK 

P(7J»(P(1)*PI3)*PI2J*PI4)J/PK 

P C ?! ) = CP (? )*P( 3 )-PU )* D C41I/PK 

°( 0 ) = l-R *( t»y*G)**2 + Ct 1 )**2)/( ( PM*G )**4+Q { l ) **4- ( PM*G*y ( 1 ) )**2 ) 

P| 10) = (-K *(PP*G)**2 + U2)**2)/( (PM*G)**4+U(2) **4- (PM*G*Q( ?>)**2 ) 

A ( 1 , 1 ) = 1 , 0 
A f 1 » 2 ) =C • 0 
A{ 1, 3 )=1.C 
A f 1 » 4 1 = C « 0 
At l * 5 ) = l . 0 
A (1,6) =0.0 
A( 1,7 1=0.0 
A(1,8)=C.G 
A(?,l)=CnSH(Q(l)*F) 

A(2»2)=5I NH U(1)*F) 

A ( 2, 3 ) =COSH ( G (2 )*F ) 

A t 2 » 4 ) = 5 1 NH < G { 2 ) *F ) 

A(2»5)=C0SH(G(31#F)*CC:S{G(4)* : F1 
A C 2 , 6 ) = S I NH t C I 3 ) *F)*CCS(Q(4)*F> 

A(2,?)=C0SH(C(3)*F)*SIN(Q{41*F) 

A ( 2 , 8 ) = S I NH tC(3)*F>*SIN(Q(4)*F) 

A( 3 i n=oc 1 )**? 

A (3, 2 1=0.0 
A(3,3>=<J[2)**2 
AC 3,4 1 = 0.0 

A ( 3 , 5 )=0 ( 3> **2-0 C 4 ) ■ . 

A {3, 6 1=0.0 
A(3, 71=0.0 
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i.i.R r r 0 T A T f- r f : T 


FORTRAN SC URGE LIST 


ft 

{ 3 

r ) = 


3) 

* 5 ! 4 

) 



( A 

1 1 = 

c (11* 

*2 

*A ( 2 

, 1 ) 



l 4 

? ] = 

0(1) * 

*?*A ( ? 

.21 


• 

( A 

1 ) = 

•At) * 


*A (2 

, :•) ) 


A 

( 4 

A 1 = 

0(2) * 

J, *> 

*A (2 

, A) 


A 

f A 

5 1 = 

A ( 3 , 5 

)* 

A(2, 

5)-H2)*A(2, 

a ) 

A 

(4 

ft 1 = 

A (3,5 

) * 

A{ V , 

ft) - i ■ ( 2 ) =£ A ( ? , 

7) 

ft 

( A 

7 1 = 

( A ( 3 , 

5) 

-■ A ( 2 

.7)+P( 2)*A(2 

,6) > 

* 

/I 

( A 

8 1 = 

{ A ? 3 , 

ft 1 

*A(? 

,fi)+P( ?)*A(2 

,5) ) 


A ( 5 , 1) = p ( ) ] 

M5, M = C.O 
a ( 5 , )=p t i'.i ) 

M 5 , A )-C.n 

M s , s i=p(y ) 

A ( 5 , ft 1 = 0. 0 

•M 0,7 1 = 0.0 

a t s, ? )=p( a l 
A ( 0 , 1 ) = 9 ( 9 1 * A ( ?. , 1 ) 

A (ft, 2 ) = tM<n*M? f 2 ) 

A (ft, 3 ) = P( 10 ) * A ( 2 , 3 l 
A (6, A > = P( 10)*A(?,4) 

A t ft , 0 1 = P ( 7 ) *A ( ? , 5 ) +P ( 5 1 *A ( ? , R ) 
A ( fc , ft ) = p < 7 1 * A { ? , 6 1 + P < 5 ) * A ( ? , 7 1 
A 1 ft , 7 } = P ( 3 ) * A ( ? r ft ) + P f ft } * A ( 2 , 7 ) 
A ( 6 » f! ) = P t 8 )*A(2,5)+P(ft J*A(2fR) 
A (7,1 ) = K( 1 ) ( 9 ) 

A ( 7,21 = 0.0 

M7»3) = G(2)**2*P(10) 

M 7 » A 1=0.0 

Af7,5)=P(7)*A(3»5)+P(21*P(5| 

M 7, ft)=C.C 
A(7,7)=C.O 

A[ 7,R)=P(R)*A(3,5)+P(2)*P(6) 
MR, 1 ) = M A » 1 )*P{ 9 ) 
A|fl,2)=A(A*?)*Pt5) 

MR,3 ) = M',»3)*P( 1C) 

MR, A )=A(4,4)*P(1C) 

A( R,5) = P(7)*A(4,5H-P(5)*A(M6) 
M8,M = Pm*A(4,6) + P(5)*A(4,7) 
A(8,7)=P(8)*A<4,6)+P(6)*A(4,7) 
A ( R,R )=P[R )*4 (4, 5 )+P(6 )*A(4,fl ) 
R( 1 , 1 >=—1-0 
B ( 2 , 1 ) =- l . G 
9 ( 3, 1 1=-(PM*G)**2 
B ( 4 , 1 ) =- t PM*G 1 **2 
0(5, 1)=0.0 
«( fe, 1 )=C.O 
B ( 7 , 1 ) =C . 0 
0(8,11=0.0. 

N = R 
: N = .i 

CET£Rr=0.0 
15 DO 21 .1 = 1, k 
2i IPIVCTI J)=0 . 

30 DC 55C 1 = 1, N 



o n r> r> n n 


SCU :r.r STATEMENT 


F OPT RAM SOURCE LIST 


r 

C 

r. 


SEARCH Fi ik PUTT FL r - T 
6 f, 

nrj if 5 .!=!,?: 

5 r I F ( I p 1 VC T ( J I - I ) 6 0 ,105, 60 

60 rn ICC K = 1 , .*•; 

70 ir-f iPivnnKi-i iso r 100, 740 

-if) n ( AilSI AM AX ) -A 05 U (J , K ) ) ) 35,100,100 
n r > iunu=j 
•90 rcn. c ‘-'= k 
oo av a <=;.( j,k> 

100 CONTINUE. 

IGF CC-MIkiJF 


C INTFRCHAUGil ROWS TO PUT PIVOT FLOHHNT ON RIDGNAL 

r 

1 in 1PIVCTI IDjLUK)»IPIVGT( ICOLlJP ) + 1 
1 10 rot IflCw-ICCLUKJ 140,260, 140 

160 CcTlPP = -DFT!:RP 
150 DO ROC L*1,N 
160 SW*P*A( IP.CV.’jL ) 

17 0 All R G V- , U =A (I COL IJ M , L ) 

2C0 AMCOLOVrL INSWAP 
705 I F I M ) 2 6 G » 2 6 0 » 7 1 0 
210 DO 25C L = 1 * Pi 

2 20 SWAP=P< I ROW , L I 

2 30 D ( IRGl. ,L)=P( I C C L U V » L 1 

2 50 Q( ICCLUK,L»«ShAP 

260 I r.DFX ( I , 1 ) = IROW 

270 IKOFX ( I ,2 UICCLUV 

310 PIVOT ( I > = A ( I COLOR , ICOLUH) 

3 2G CETFRR = DETFRP*PIVCTU ) 


Cl VICE PIVOT ROW BY PIVOT ELEMENT 


330 A( ICCLUV, 1CCLUVU1. 

360 DO 3SC L = 1 , N 

3 50 A{ ICCLUM.L J =A I I CGLUK , L ) / P IVOT ( I ) 
355 IF (tf) 300,380, 360 
360 DC 37C L»1,M 

370 HI ICCLUV, L)=HI ICCLUM.L 1 /PIVOT(I) 

REDUCE MON-PIVOT ROWS 

3e0 DC 55C Ll=l , N 
390 IF(L1-IC0LUM)6CC,55G,600 
600 TA=A 1 LI , ICOLUM) 

620 A ( L 1 , ICOLUM ) =0. 

430 DO 65 G L=1 ,N 

6 50 A(Uf LJ=A(H,U-A (ICGLUM,L1*TA 
655 IFtM 1550,550, 660 
660 DG 50 C L=1,M 

5G0 B ( L 1 , LI =B ( LI , L l-B ( ICOLUM , L)*TA 
550 CGNTINUF 



55 




FC«TR\*J Sf'L'RCH LIST 


I 


CfiLU 


1 1 


C 

J 


6C0 nr 71 C I = L , 
oio L= f ; + i-l 

620 [F( P.TrlXlL.D-IWCEXd,?) 1630,710,630 

n30 jnnt.= ioi.:e* u,n 

64 c JCCIJ. ,^-l*:|;! r :x(L,2) 

650 f : f : 7C‘i K= 1 , ,'l 
6 60 SwA1>*MK» J??CJW) 

6 7 n MK, J R f» 1 -J ) - A I K , J c a L U M ) 

76 c uk, . jrnu;r- 1 =sv,ap 

7 0:', f.U.'O'T I !\ ! ; F 
710 CiiKTIM.ir 
7 AO SC?a I M;r 

!*K I i ! T 02 0 I (M I ,1 I , 1 = 1,*} 

S?C FORMAT! 1X,P:E15.4) 

Y = 0. 0 

TOCO T ( ll=COSH(ij(3)*F*Y) 

T(21 = SINH (U(31*F*Y) 

T(‘*)=cns ( 0 ( A ) * F * Y ) 

T(4)»5IN('J<4)#F*YJ 
t { o ) =rnsH( r j ( l ) * f * y ) 

T!6)=SIKH (CI(l)*F*Y) 

T( 7 )=CGSH(g(2)*F*Y) 

T ( F ) = 5 1 MM UM2)*F#Y) 

1, 1)*TIS) 

V. ( 2 ) -= M ( 2 » 1 ) * T ( 6 ) 

W(3)=f;n f 1)*T(7) 

W{4)=LM A, l)*T(fl) 

W ( 5 ) = IK ( 0 , l ) # T f. 1 ) # T l 3 ) 

W ( 6 ) = <5 ( 6 « 1 ) *T { 2 ) *T ( 3 ) 
h(7)*R(7tn*T(l)*T(A) 

V< ( B ) = P ( >3 « 1 ) *T ( 2 ) # T ( 4 ) 

K < 9 > ( 1 ) +W I 2 ) +k ( 3 > + U ( 4 ) +W ( 5 ) +K { 6 1 +W (7 J + W ( 8 » + 1 . 0 
WW = 5 . 2* Vi I 9 ) / ( PH*G ) **3 
. XV(i)«-hC9) 

XfM 2 ) =S ( 1 ) +*2* ( R ( l * 1 ) *T < 5 ) +R ( 2 1 1 ) *T ( 6 I ) 

X* (3)=U(2)**2*(B{3,1)*T(7)+B(4,1)*T(8) ) 

WK=rt ( 5 1 +W l 6 ) +W I 7 ) +!« I 8 ) 

X M ( 4 ) -( y ( 3 ) **?-G ( 4 ) **?. )*WK 

XM(5>s2.*<j(3)*G(4)*(-P (5, L)*T12)*T(41-B(6t 1 |*T< 1)*T(4)+B( 7, 1)*T(2) 
?*r(3i+B(6fii*T(n*T-{3n 
ZM= (Xl^(2 )+Xf'( 3 J+XPM* J + XMC5) 1 / ( Pf'i’f'G ) 

(XP ( 1 |+. 3*ZR+1 .0 )/ { 2. 1*PM*G) 

YMM=-(ZR+.3*XM(l)+i. )/ I2.l*PM*G) 

X YP { 1 )=Q (ll*(B(l»l)#T(6)+B(2ii)*T(5) ) 

XYV { 2 1=0 ( 2 1 * f Bl 3 » 1 1 *T 1 8 ) +BI 4 » 1 ) *T l 7 ) ) 

X Y I* { 3 V=Q I 3 } *( E { 5 , 1 ) *T ( 2 ) *T ( 3 ) +& ( 6 , 1) *T C l ) *TI 3 ) J 
XYl»14) = Ct31*< B(7»1)*T(2)*T(4)+B(8, 1 ) *T { L ) *T ( 4 ) ) 
XYK(5)«C(41*(-0(5.Vl*Ttl)*T(41-B(6»n*T(2'l*TC4)l'- 
XYM ( 5 ) = C C 4 ) * I-B ( 5 » 1 ) *T ( i ) *T (4 )-B ( 6. » 1 )'*T ( 2 1 *T ( 4 ) ) 

XYM ( 6 ) =Q{4 ) * ( B { 7 1 1 ) *T ( l ) *T ( 3 ) +B (8 » 1 ) *T (2 )*T( 3 J > 

XYM(7)=iXYH{l) + XYM(2) + XyM(31+XYM(4l+XYH(51+XYM(6j 
, XKM*-XYMt7l /(3iO*(PH*G)**2) 



FORTRAN SOURCE LIST 
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SOURCE STAUSFEfcT 

? 1MI I T /■ 0 C 0 , Y , V. v, , X F‘ V , Y f-' , X K P 

* 2CCC RiRr'.vn IX ,^f:l c i.4 ) 

Y = Y + . 1 

3 1 F ( Y - l . 0 ) 1 4 0 Q » Vi C 0 i 2 1 C 0 

S ? ICO F=F+ . ? 

7 l F ( F- 2 . )73*7fi,2?CC 

;) 22 CO X*X+5.0 

1 IFJX-50.0 >2,2,6 

? 6 FM=pr+?.Q 

3 lFU ) f'-r.J'»*AiO 

4 <5 5TCP 

5 N n 

ITWAR ASSEMBLY 


ESSAGFS FOR AMOVE ASGFNKLY 


o r- 



0 » IHV-TC CL'Mf; t'CRECK 

1 SlIll^rLTiriEi CLfMC{A|XRiXI ) 

? 0 [ '■ f: N I G? A ( A ) , X R ( 3 ) I '■ r3(3) 

3 I PATH-? 

A FX=l./3. 

- IQftn F0«r:AT(5l'l^.2l 

ft [ r ( a ( M ) 1Q06, IQOA , ICQft 

1 100 A XR (11=0. 

0 GO TC 10 3/, 

1 100ft A? = AI 1)**2 

? ■;= \ 7.1 .'"A?-* A I A )-9.*AI 1 ) *A U)*A l 3)+2.*A l? )**3 ) / I 5/, * *A2*Al 1 ) ) 

3 I F ( 0 > 1010 » 1000,1014 

A 100? 1 = 0. 

ft GO TC 1032 

ft 1010 

7 I P A T h = 1 

0 10 1 A P=(3.+A(I.)*A(3)-A{2)**?)/(9.*A2) 

1 A “■{ G = P X !5 l! 3 + * 2 

? IF (ARC) 1016, 1013*1020 

3 1016 Z=-2.*S0liT(-P)*CC5IATAN<SyRT(-ARG)/ti)/3.> 

A GO Tf 102 8 

5 10 1H Z--2«'H:**FX 

GC TC 10? a 

1020 SARGsSORT (ARC) 

0 IF (P ) 102? , L02A, 1026 

1 102? Z=-(C+5ARG)**LX-(C-SAHG)**EX 

? GC TC L 02 3 

3 10 2 A Z=-(2.*0)**FX 

A Gu TC 102b 

ft 1026 Z»(SA»«-(n**EX-< SARG*U**EX 

6 102 3 GC TC (1030, 1032) , I PATH 

7 1 0 3 C 1 -- 1 . 

0 103? XKI 1 ) = ( 3. *AU )*Z-A(2 ) I/I 3.#A{ 11) 

1 103A AQ(1J*A(1) 

2 All (2)=AI? )+XR(l )*A(1) 

3 AO ( 3 ) ( 3 ) +XR I 1 ) *AQ I 2 ) 

A CALL CUAD(A0>XR(2) ,XR(3) , XI J 

ft return 

6 END 

IBMAP ASSEMBLY CUBIC 


FSSAGES FCJR ABOVE ASSEMBLY 
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FORT TAN ", i-'Cr i f'-r 

1 SOURCE STATEMENT 

) $ I B F T C OUAD fiODECK 

SUBROUTINE OIJ4U(A»XRl ,XR?, <1 ) 
j DIMENSION A (3) 

\ ?.h FORr'ATf 5F14.3) 

y X 1 =■" A ( 2 ) / ( 2 • * .4 ( 1 J ) 

> DlSC=Xl**2-AI3»/Al l) 

, IF(DISC)?7,2«,28 

r 27 X 2 = SQRT (-DISC) 

) XR1=XL 

XR2 = XI 

! X I = X2 

i GO TO 29 

y 28 X2=$QRT (DISC) 

i xru=xi+x? 

, XR2=X1-X2 

' X I =0 . 

> 29 RETURN 

END 

I PM AP ASSEMBLY iJUA r ; 


rSSAGES FOR ABOVE ASSEMBLY 

IBLDR — JO 8 0000 DC 


IS BEING ENTERED INTO STORAGE. 


TWO iiUvii SIMPLY SUPPORTED <* TWO iiuG^S CLr-MP^u 


F i ; F f" 0 A ‘ , V UP Cl LIST 

:• SO'l/XCt'. 5T/.T'::r»lT 

I fi h T 0 ..OOFCK 

1 0(b) ,-T(4) , Ar;(4) , ■- 0 { 4 ) , i I ( 4 } ,A0UJ , XP(4 ) 

;• u i ■■'-■ < r > I cm r h :• ta { 4 )^?K),oi{aj,f- c 4 , , f 1 1 a , ; ) ,p r ( « ) 

5 1)1 .‘1 •.: '-J ii 1 0' i I PI 7' I T ( :j ) , A ( O , ) * ; { 3 , 1 ) , I MO V ( b » ?. ) . I V 0 T ( H ) 

4 p ( i ) , w(9J , i 01 r x*<(b ),o(Mi xvMm 

C F - / A 

5 ! - y -t 

6, r>K Ii-T'-'. 

7 * F' A T ( I H 1 ) 

0 4 F. ~ -J <: 1" 

1 ?‘=0„C 

2 2 r- = I - c 

3 7 1 = 7 = 1*16 

l FACTF!=Pir*PM 

5 7 = 0 > 3 

6 U=i te *{ i« 0-V**2) *X**2 

7 c ( i ) = i « :? 

0 C ( 2 ) =~4o 0* ( FACTR )**2 

1 0 ( 2 ) =<:* 0* ( FACTR )**4+L 

2 C(4 J =-2« *( FACTR ) **2* ( ( 2* *FA'CTR ) **4+lJ*P, ) 

3 C ( 5 ) - {FAC TR )'**4*( (FACTR )*«4 + R'**2«U) 

4 A3=C(2)/CU) 

5 A2=cm/C(i) 

6 Ai=Cl4)/CU) 

7 A0*C(5)/C(1) 

0 AT=A3 /2o 

1 AC(i)*l. 

2 AC { 2 )=-A2 

3 AC ( 3>»Al*A3-4.*A0 

^ AC (4) =A0* (4®*A;2-A3**2 )-Al**2 

5 CALL CUBIC UCt-RT'.RTU 

> I F ( R T I ) 2 0 > 1 0 r 20 

7 10 I F ( RT (1 )-RT ( 2 ) ) 1 1 1 12 » 12 * 

3. I I KT ( 1 ) -RT ( ) 

L 12 , IF CRT (1 HKT < 3 ) ) 13>20»2Q 

■ . .o, rtu 

5 ' 20 b-RT ( U/2® O'- 

t . 1 F ( b^*2r-AD ) 2? r 22 >24 • 

c 22 ' ; 6 ?»Q. ?. \ . 

rV ;;••• qo' Td :; 25 '- ; j;o v:;' y V ; ..v '• 

M v ' •: ■' . •••.' '.'.'.O' 

C A=~ ( «?i/2 e-A r*o ) / 1> 

l 25 ALi(l)-I« 

) A& (2.) ®AT-CA . 

» AH(3)*0-0 ; , f . „ . . 

f... CALL CUADCAQ.t XR( 1 ) t 2 ) » >0 ( 1 ) ) 

> . •’ s?QU 

! :• bQ(2)-AT+CA r ; ' 

.call; ' cOao- 4 CLttt'Xijtf 3> vXRtA 

«;■. 'v.\- ' t £ Vy-./- ■: . . 

c:*ov- f xr it 4 ¥ U$V : . ‘ '■ ■ • 

,P1UHT?T;*X ' V-- • ;•"• • •• ' ; 
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t ?,pf 

7 7 FG! : .‘-A f ( 1 X i P 7 « .? ) 

fj 00 i 7-1,4 

1 IJ/xsin) 

2 TH': f A ! T ) =4 r Mi C P P f I ) ) 

j P 1 = 4 , # m T A ! i ( 7 - ) 

? * ’ : < r, ( i J =S(J«r ( Xh( [ ) **2 + XI ( I )*«2) 

5 i'. KI) = r,CKT ( V>. P ( I ) ) 

FR ( I , 1) = RH I I ) *COS ( PI+C. t »*TH;;T4 ( I ) ) 

7 F I ( I , !) =!<>’; ( I ) »SI N ( f 3 1 + G „ j*THi: T 4. (II ) 

r ! F I( I , 7 ) = <h ( I ) *f,IN( *P I+r ; , 5*ril r ; 7 A( I) ) 

^ * F 1 ’ ( I i 2)=XU( I )*CQS(2«* 1 .'’ I + 0*5*Til !? TA( I ) ) 

’• PkH-jrfi, UTRd, j) ,I=i,4),J=l,2) 

^ P! '- I f 3Tiii((FI(I l j) t l = i t ^j t j~]_ J 2) 

5 a F U R A T ( 1 X , a C 1 4 * 4 ) 

C INITIALIZATION 


6 

7 


Q(! )=Ff< (1,2) 
C( 2 )=FX( 2 , 2 ) 
C ( 7 ) = FX (4,2) 


1 Q(4)=F 7(4,2) 

2 G =3*1416 

3 Pf 1 )'="■* 5* ( PM*G )’**2+ Q(3)*^2- Q(4)»-s2 

f i P{2)=2**Q(3)*Q(4) 

5 P ( 3 ) = (PH*G)**4+Q( 3) «*4 + Q ( 4 ) **4-6* » ( Q( 3 

2* ( PM*G)**2 


u'-'C; l i x r 


)*Q( 4 ) ) ** 2 - ( 0 ( 3 ) * * 2 —Q ( 4 ) * * 2 ) 


^ P (4 ) -4a *Q ( 3 ) **3»Q ( 4 ) ~4»*Q( 4 ) 3*Q ( 3 )-*2. * 0 ( 4 ) »Q ( 3 ) ■» ( } **2 

7 PK = P ( 3 1 **2+P ( 4 ) **2 

C- P'(5) = (P(4)*P(1)-P(3)*-P(2)')/PK 

1 P(fc)=(P(2)«P(4)+P(l)*P(3))/PK 

l P(7>=(P(1)*P(3)+P(2)*P(41 >/PK 

3 P( 8 ) = (P (2 ) *P (3 )-P (1)*P(4))/PK 

5; ,P( 9) = (-«.5* l-P.M*G)*.*2+Q /( ( PW*G)**4+Q( l)**4-( PM*G*Q( I) )*»2)' 

i> P ( 1 G ) = (-* 5* ( PM*G ) **2 + Q(H ) ®*2) / ( ( PM* G) *®4+C/ ( 2 ) »*4- ( PM*G*Q ( 2 . ) )**2 ) 

& ■ ■ A(1,1) = U0 

? M 1 , 2 ) fOi 0 

9.- A 6 .X ' :■ •• 

i ■'■■■ A'(i ; ..,'4 : ')FO»0'. ■ 

A.:(TV5.)=1^0,:;: •• 'V-' 

* / ■ ' 7. ./ •••,.■ • . . • . ■ . . , ■ vP 

f-‘>V7-.;v, 7- y.i •••••• 7.^7/; : v; .. ./ .7..' 

f Mi’ 4 f= 0 , 0 '*' 1, '' ' 

S A(2>1)bCQ5H(Q(1)hF) 

i ;a<2,2)=sinh (ou)»f) 

> M2»3>*CCISH(QU)«F) 

t AUV4 )=5I7?H l Q ( Z ) *F ) 

? . M£*S)«CQSH((i(3)*F}.*CCS(G(4)*F) 

* M2*6)b&1NH {G|.13)*F)*C0SIQ(4)*P) 

A< 2>.71«CaSrtfQ-C3 )V-ri 4 MF 1 

V.:- 

■; :■ v .. 

£ . ' ... ••:••• : 

f- if /--- 1 •• 

?«.• •: • - 



O ^ -V* & ^ O -£> 0* l\j »-• Or-J 0> 


M 5 » 7 1 = * 

:( 

^ ) 

M ■ , i: ) = ' 

■ s 

r> 

/■ ( 1 i =c 

{ 

i)*A(? ,2) 

/- ( ^ i 2 ) - 1 

,( 

1 ) * A (2,1) 

A f A , ■?. ) =< 

. ( 

) * m , a ) 

A ( * , A 1 =C 

,1 

2 ) «• A 12,3) 


M 4 , 

A 14, 


-!.{ 3}*AI2 , A)-Q( 


) * & ( 2 , 7 } 
) * A l 2 * M I 


A ( 4 »M«0( 3 )*A(2*5)-C( 4 ) *A 1 2 , M 1 
A « ( 4 , 7 ] =(; ( 3 ) *M 2 , >-; )+(■ i I, ) * A( 2 , 5 ) 
A ( A , f! ) =g ( 3 ) * A ( 2 f 7 1 +G ( '•> ) * A { 2 , M 
A(t,l)=P<9) 

A f " » 2 )=C,0 
Ml ,3) = P(1C) 

M 5, A J=C*C. 

A15, 5)=P{7> 

A ( 5 , 6 ) = G» 0 
Alt *7)»C,0 
A(S,P. )=P<8) 

A(t,l)=P(9)*A(2,l) 

Alfi,2)*Pl9)*At2,2) 

A(8,3)=P{10)*A(2,3) 

A(6»4)=P< 10 ) *A ( 2 , 4 ) 

A16»5)=P{ 7)*A(2, 5J+PI 5)*A(2 ,8) 
A ( 6 , 6 ) = P ( 7 ) *A ( 2 , 6 ) +P < 5)* A ( 2 , 7 ) 
A(6»7)=P<8)*A(2*6>+P<6)*A(2|7) 
A ( 6 i 8 )=P ( 8 1 * A ( 2 , 5 ) +P < 6) * A ( 2 , 8 ) 
P(ll)=0l3)**2-Q(A)**2 
P(12)=QU)**2«A(2il) 

P( 13)=QU)**2*A(2»2) 

■PI 1A>*0(21**2*AC2-i3) 
f(151 =Q12')**2*A(2,AJ 
P(16)=P(ll>*A{2,5)-P{2)»A(2,8) 
P l 17 ) =P,rii J.»'A ( 2, 6 HP ( 2 j*A C2 , 7 I 
P ( I87“P( 1 1 ) *A ( 2 * 7 ) +P ( 2 )* A ( 2 » 6) 
P I 19} •=? (,1 1 1 * A l 2 » 8 ) +P ( .2 1 * A ( 2 , $ ) 

■ . A ( 7 f l l=G(l)*fr2*P(q} . 

A ( 7 ,37=:0{ 2 ) * *2>P ( i.0 ) 

+p(2>*p(5> 

v ■w**)?p0® . iv •• 

A ( 7 , 8 ) = P ( e ) #p ( U ) +P ( 2 1 *P l 6 1 
AUtl)*P(12) *P<9) 

M6r2)=P(l:n «P(9) 

A{ 8,1 }*P ( 14 ) *PU0) 

• Af3*4)-P{m *P<1 CJ • 

- - AC fi.t )*>( 16). *P m*P( 19.) 

AMr»&'**Rf 7I*Pl IT ) +P ( fJM.Pl 18.) 

- ■ . i ■ A(M1*MB):»Pn7) ->fr{fc>*PI i&v 

" 1 6f +P ( &)*Pf 19) 

• Mi } ■$ *. •• /*. : • •• 

• r'-v- • .. ": 



P-.’-P r 0 LI IT 


6S 


F.n,MCf 



H< I ) = f! = v 

L 

0 ( t , ) = C,0 


('(*,.'.} = <i> ij 

t 

, M ■■' t j. ) = ■- B i- 1 

■, 7 


'0 

f / = * 

rl 

f J c I ./•/" = '.,« ‘J 

rz 

Of' <:■ J=l,^ 

n 

1 IPIVCTC J)='J 

rj 

V; r.;L 1 = 1, 


c 

C f»-_A3«TF mk PIVOT r.L zv '.: 1 l T 

c 


? 6 

✓ r* 

*T 

A I' / - X ~ C r, 

1 7 

■’T 

[jU j. ..* .j . ! — j. t i 

A 

50 

IT { IPIVOTi JJ-1)60,105,60 

51 

/ r- 

o u 

DC ICC K=1,N 

•. *•* 

70 

If ( I PIVCT ( K ) -3 380,100,740 

;3 

y o 

IF (ARS(AMAX)-AHS(A( J,K) ) ) 

54 

85 

IRDW= J 

55 

. 90 

ICULUP-=K- 

56 

95 

AtfAX=A{ J,K) 

57 

IOC 

CONTINUE 


1 105 CONTINUE: 

c 

C INTERCHANGE' ROWS TO PUT PIVOT ELEMENT ON DIOGNAL 

C 

,3 no IPIVOTI ICOLUM) = IPIVOT nCOLUM) + l 
.4 130 IF ( IR.CW-ICOLUM >140,2.60,140 

.5 . 140 DETERS- DETER* 

6 150 DO 2 G C L»1,.N ' - 

.7 . 3.60 SWAP=A( IRON, L ) 
fG 170 AI IRDK, L ) =A( ICQLUHrL > 

'll 200. At ICCLUM^L INSWAP ,' : - 
%■ ■ .205. IFT*) 260,260 ,210 ' : . 

A- ■, -210 \Dtl.v25C "L'-i . 

[ 5 .' 22-6 

IB 230 R(7T$C^U')=Bi(tiCWR,tl' ■ 

;t -.v:; o- 

ix ■. : :■ . : • : ’ - : 

ty, . ■ 2 . 1 T 0 - ••/»•. .• :. ' • . iuV .■ .?:• 
'p.ivd.ui»A i iconw i-:v .•••;:• -j.- •. <: 

c 

C DIVIDE PIVOT *GW BY PIVOT ELEMENT 

C ■ ■' 

15 3.3C A( 1CC:LUK> ICHL'JPi J “i* 

1ft '340 »C 35C- t>U* 

\T 350 A f I CCtOM.rt ) -A t r COL UP- , L ) /P I VC t U i 

il. ... 383 :;IFt^me*3«0,^6Ti • • ' 

\% '' ' 36 R' .DC.- 37 £kL-*ii&’ ■■• 

f> _ : \&70 '••« U fP I V GT l -t I 

t d oucie •• fsf'ofe— vqx . •: : /•“ 
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f=? jp r r .v-i ecu- cl ltj 


D'JV cr s T A T i- 7 i ; T 


0 7 -jC L.T. 6. :: LI =!»(•: 

1 6 - v 0 I F ( L i - I C 0 L 0 ) L 0 , 5 5 f; , A A Q 

i7 4C-J T/.--MU , ICC LU ?f ) 

t-l ( : A (LI » ICC LUf)=0i 

;j. av iv. l = :»m 

: 2 -7 6 r .'. -Ml] 7 L ) = b ( L I , I. ) - A {] C 0 l AM , U * T •*« 

i 4 f * L 1 1' ( 7 ) 1- T> 0 i 5 : T . 7 > A £,. 7 

• 5 fc a fj f;fj CCC L = 1 » V 

■A. 5CC HtL; » L ) = i ( L 1 j L ) — C ( I C L L IJ?' , L ) * T A 

: n 5-ji' C • ’i-:T r-Uti 

c 

C I uT'-P CHANGE COLUMNS 

c 

»3 600 DC 717 

V* 6 1 <: L-N+l-I 

>5 620 IF t INCEX(l,l)-INDEX(L,2) ) 6 30 , 7 10 , 630 

;fc 6 30 JkGW= INDE X ( L » I ) 

■7 600 J C C L U M = I N D E X ( L ,2 ) 

10 600 DO 70S K=1,M 

fl 660 SWAP»A(K, JROU) 

'2 670 A ( K j JIUJW ) = A ( K » JCOLUM ) 

S3 700 A (KfJCCLUM )=SWAP 

fA 705 CONTINUE 

rfc 7 1 C CONTINUE 

;c 7 AO CONTINUE 

)1 PR I NT &20 » ( R ( I » 1 ) j 1 = 1 » 'ii ) 

lb 82 0 FOP. MAT ( 1 X » 8E'l 5* 4 )' 

)1 v = c»c 

.0 19GC T(1)*CQ$HU|C3)*F*Y} 

11 T ( 2 ) = 5 1 N H (0(3)-* F*Y J 

.2 T ( 3 ) =COS (0(4 )*f *Y) 

.3 Tt4)=SIN(Q(4)*F*Y) 

A T (5 ) "CCSKIQIl }*F*Y1 

:5" t{«)-SlNH (OU)*F*Y| . 

.6 T ( 7 ) = C(JSH (012) *F*Y) 

.7 . T t © h= S:i to ( Q (2 ) 

so . r 

fl . W12> = fi(2il)#T(6) 

\z : ;v;i3:) = B(3rn*T(7) 

!t> 'VfaisttfsvwT (2 i'i'ffsi’ •-•••-••••' 

:6 0 f. n^e(7,i)«-m)*HA). 

17 . W t ':• ) = H ( 3 » i ) «T C 2 ) *T ( A ) 

10 U{0)= : ka)+l«C2)4W(3]+W(4)+W(5)+W(6i)*W(7)+W.I8) + l»0 

is npni^i?) ••• •' 

I A , . ' ■ XM t i ) *Q { (313 1 1 1 * T ( 7 1 + & C* » 4 ) *T 18)1 
|5 ;l g;l+« H-) +.Wf.n f»'t $ f >' 

ii>- (4 } * « 2 1*WK.; 

it . :& ( tn *■?:* *9i -U *« i a ) n **•« t $ V i > *T:( 2 > *t.i 4 ) »t ( a j +b ( 7 . »-x ) * 

c 1 . *. .t* $ if 4 t * €b t / < • 
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r nu':C" GT 7 - 


::T 


%• t j - : C i 


[ i. 

L ’J 

14 

L 3 

it 

,7 

it 

?1 

: *7 
_ 

?4 

?5 

56 

n 


, I [J j- r c C 1 ' *? It ! j 0 o ; ; C H 

LUG'- C L T I ’ S l CL tic (A, /?•:,/ 1) 

U ! A f 4 ) ,XMi) ,iC(? ) 

I ■- u I r ■» l 

tx = l « 

it L0 FC-f'J-'iiTfiFiA* 2) 

IF (4 (4) ) 10 L 6,1004,100 4 
1 C 4 X!'; ( I )=G« 

GC it. K34 
IOC Ei 4, : ' = /-(;.)* **. 

0= (27,*62« MAl-'J.m I )*A(1 ) .A I 3 ) + 1 R * { 2 ) * 
IF ( 0 1 U> i. 0 i IO'jC , 1014 

1008 i - 0 r 

CC TC 1032 ■ 

icic os-o 

T P A T H = 1 

l C 1 4 p= { ?, « *A m * A ( ? I - A ( £ ) * *2 } / ( 0 , *A 2 1 
A ft G = ?*»•. + 0 » *2 
IF(ARG)l0i6, 1018*1020 

If;!/:. Z=-2*»SL'RT f - P ) * C 0 S ( A T l. N ( SQRT(-ARO)/Q) /3. I 
GO TC 1028 
10 If Z=-2« *C**FX 
GO TO 1028 . 

10 2 C $ARG= SORT ( ARC ) 


10 I F ( P } 1022 1 1 024 » 1 026 

51 1022 Z's-(C + 5AKG)«*EX-(Q-5ARG)«»EX 


52 GO TC 102b 

13 1024 Z=-{2 e *0)»*SX 

14 GO TC 10.28 

55 1026 7.= ( 5ARG~G)**£X-(SARG + 0)**EX 

56 1028 GO : TC (1030,1032 ), IPATH 


57 103 O' Z=-Z 

jO 1032 XR(1) = (3**A(1) «2—A{ 21 ) / ( 3* *A( 1)1 
)1 1034 AG( 1 ) -A(l ) 

(2 A6 ( 2 1 = A (2 ) +.X R ( 1 J * A { 1 1 

f3 ' AO l 3 ) p:A ( 3 ) +XRI 1 ).*A;Q 12) 

i4 CALL OUAD ( AO, XR< 2 1 , XR (3 j , XI ) 

».5 " • • RETURN 

•6 ■ END. • ' • 


I SWAP ASSEMBLY 
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F'-jH o..v; S* ' C'. 


r_ " ' 

Li r, .• LTA T ; v - 

:.t 

r ; 

^ j fj -TO loo ! -;r;r» ^ c 


T 

L 

jL’-'i- OLT J.'iE OLO- 

U ( A t X h I 

? 

o i i - s f ii;:-i a n ) 


- 




■< 1 =- A ( ? ) / ( , * A 

(III 


s; I iC = x { 

i)/4(i) 

r 

IF (0 I EG )■■: 7, 2.i, 

2 n 

7 

2 ( XZ*Vm r < — 01 SC 3 


10 

A *' i ■“ / SL 


u 

;< ? : = X 1 


3 >3 

x t = x ::• 


1.3 

! 3 :j t r: ?7 


14 

/■' / 2 ~ S C, )•; T ( u i G ^ ) 


IE- 

X?-. l=XI + :<2 


16 

aV; ; = xi— x 2 


17 

/ x ~ O * 


20 

2 9 BETUHiN. 


21 

!:ND 



SAGES FOR ABOVE ASSEMBLY 


I3MAP ASSEMBLY 


IBLDf! — JOE 


¥ 15 BEING ENTERED INTO STORAGE 


LI T 


OU&D 


00 GOO 0 




















10 



V«*2.4 

V r ®* s 0.8 



Y/b 

io 3 ^^ 2 

K&V^iir 


10 3 Msy/gJ.T& 2 

0*0 

0.00 

-10.61 

0.00 

-10.18 

0.1 

IS. 66 

-18.16 

—21,66 

-15. 10 

0.8 

is.es 

-18.86 

-81.13 

5. 4? 

0.3 

1,03 

-16.06 

- 16.02 

16.02 

0.4 

* l.BS 

-14.14 

-12. SI 

11.08 

0.6 

- 6*91 

-13.43 

-11.08 

0.00 

0.6 

— 1.06 

-14.14 

-12.18 

-11.08 

0.? 

1.83 

-16.06 

-16.02 

-16.02 

0.® 

16.86 

- IS. £0 

— 01, 13 

- 6.4? 

0.0 

1S.66 

-18.18 

-21.66 

26.10 

1.0 

0.00 

-10.61 

0.00 

10.18 

^piiFj»i0 

fe/«*i,6 

*/»K>,0 

R l/Bj^O.8 


¥/b 

io 3 tfVXr* s 

10 r «j/£itir 


lO^y/gm 2 

0.0 

0.00 

-10.61 

0.00 

- 6.16 

0.1 

18.31 

-18.30 

-cr.so 

- 9.94 

0.2 

10.46 

-11.09 

-20.0? 

13.11 

0.3 

* 3.94 

-14.61 

-14.61 

20.09 

0.4 

-11.12 

-18,60 

-11.49 

18,98 

0.6 

*23.00 

-11.90 

-10.60 

0.00 

0.6 

*11.12 

-18.60 

-11.49 

— 18.98 

0.1 

* 3.94 

-14.6? 

-14.61 

-00.09 

0.6 

10.46 

-11.49 

-80.0? 

-13,01 

0.9 

18.3! 

— 18# 38 

-22,80 

9.94 

1.0 

0.00 

-10.61 

0,00 

6,18 


y®»i*8 

*/«**0.e 

H a/V°-* 


Y/b 




10%*y/iU11i^ 
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« e /W*l*"K> 

Y/& 

0*0 

M 
0 . £ 

0*3 

0*4 

0.5 

0*3 

0*7 

0.0 

0.0 

1.0 


X/b 

0.0 

0.1 

0,£ 

0.3 

0.41 

0.8 

0.6 

O.jp 

0,0 

0,9 

1.0 


V* 

0.0 

oa 

0,8 

3 

4 

5 

6 
7 
0 
0 
0 


b/®*s,Q 

*/ a *0,6 

R|/Rg®0.5 


ICptih/Aiw 

io*Wg«h s 

30^%/ sy - lb ^ 

10%* y / 

0*00 

-10,61 

0,00 

- 0,73 

10.66 

*18.3? 

-82.90 

• e ,3£ 

-S *6? 

-15.59 

• 17.33 

S&.S9 

86.66 

-18,19 

* 18.06 

20. Cl 

43.86 

—10«6 S 

- 11. 3l 

11.46 

49. 88 

• 9,93 

-11, £6 

0.00 

43.66 

*10.48 

-11.31 

-11.48 

86.66 

*18.19 

-18.06 

• SS.81 

8.63 

-15.69 

-17.33 

-Sir, 89 

10.6 V 

-10.31 

•22,95 

0.30 

0.00 

• 10*6 1 

0,00 

0,82 

ty' 6*$.0 

*/*v*0. 8 

sy.^0.5 


10^»V/l6 2 

10%V^7to® 


10%* y / 


0.00 

- 10*61 

0,00 

3*11 

—14# 3y 

-16*76 

0.84 

*16.83 

—19, 14 

4,8; 

-25,99 

-16*63 

0, 31 

• 24.97 

— l3*89 

J.iO 

«*24i,50 

— 10,00 

0.32 

-10.97 

-13.89 

4.92 

-16.99 

• 16.63 

0,84 

-16.63 

-19*24 

3* 11 

-10,67 

*13,76 

0*00 

-10.01 

0,00 

fcj/ 4$®1» 5* 

*/«*£• 3 

VK**e 

XQpt&j-rtP 


® JO%/l 


0 

0 

0 

0 

0 

0 

0 


(J-.CO 
0*86 
» 1. 49 

• 0.76 

*16.06 

• 10.04 

*16.00 

• 6.00 
• 1.49 

0.46 

0.00 


s.«e 

а. 0i 

1.39 

4.61 

4.44 

0.00 

4,43 

4.61 

1.59 

0.31 

б. W 




io fi Mnr/snn s 


-20.61 

0.00 

7.63 

-26.23 

■ — IS. 76 

- 1*83 

-26* 21 

-28*74 

ft,9S£ 

-14.39 

-14,89 

4.1*90 

-2C.00 

-11,79 

9*60 

— IS, 06 

-10,71 

o.oo 

-10,80 

-11.79 

- 0*©0 

-14*36 

- 14*69 

-n.90 

-26*21 

-38.74 

- 8.92 

—16*13 

•16,06 

1.86 

• 20*61 

0*00 

- 7.80 
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V®*1,4 

s/e*0.8 

?y s v*o.s 

y/fe 



10 r %/e/Th 

0.0 

0.00 

-10.61 


0.1 

1,66 

- 16 . ri 

0,00 

O.f 

0.74 

-18. £6 

-19, n 

0.3 — 

£1.83 

-ir.ee 

-17.91 

0.4 

33.04 

-11. 16 

-13.64 

0.5 

37. £0 

-10, 67 

-10.06 

0.6 - 

33,04 

-11.16 

-10.06 

0, 7 - 

£1.83 

-ir.ee 

-10.06 

0.® 

8.76 

-15.S6 

-13,64 

0.9 

1,66 

— 16 « £4 

-17.91 

1,0 

0,90 

-10.61 

-19, 24 

0.00 



V®*0*© 

?l/P£*0.5 

V/fe 

Itshb/JLivr 


ao r %/^'m r 

0.0 

0,00 

-10*61 

0.00 

0. 1 

- 3.98 

-10.21 

-19.76 

0* £- 

-16.31: 

-M.30 

-16.7© 

0*3 

-34.34 

-11*03 

-If. 47 

0.4 

-48.01 

- 9,84 

-10.60 

0,0 

-62.89 

- 9,3f 

-10,21 

0.6 

-40,01 

- 9.84 

-10,60 

0.7 

-34.33 

-11.83 

-12.47 

o.B 

-16.32 

-14,30 

-16,78 

0.9 

- 3.9© 

-10, £1 

-19,76 

1.0 

O.C* 

-10,61 

0.00 




11 . 

►08 

3 , 

£9 

16 , 


19 

,18 

11 

06 

0 

00 

-11 

,06 

-19 

IB 

-16 

30 

- 3 

,C 9 

-11 

06 


KAlxj?/c‘m r 

13,30 

7.752 

£1,61 

fl.40 

12.09 

0.90 

*ir.c$ 

- 21.40 
-£1.60 
- 7,71 
-13,38 



y/0 


0.0 

0.00 

0.1 

- 6.64 
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